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Abstract 

We explore electron transport properties in honeycomb lattice ribbons with zigzag edges coupled to two 
semi-infinite one-dimensional metallic electrodes. The calculations are based on the tight-binding model 
and the Green's function method, which numerically compute the conductance-energy and current-voltage 
characteristics as functions of the lengths and widths of the ribbons. Our numerical results predict that 
for such a ribbon an energy gap always appears in the conductance spectrum across the energy E = 0. 
With the increase of the size of the ribbon, the gap gradually decreases but it never vanishes. This 
clearly manifests that a honeycomb lattice ribbon with zigzag edges always exhibits the semiconducting 
behavior, and it becomes much more clearly visible from our presented current-voltage characteristics. 
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1 Introduction 

The electronic transport in graphene nanoribbons 
has opened up new challenges in nanoelectronics. 
A graphene nanoribbon (GNR) is a monolayer of 
carbon atoms arranged in a two-dimensional hon- 
eycomb lattice structure 1-4 and can be regarded 
as the basic building blocks for graphitic materials. 
Graphene based materials have potential applica- 
tions in several branch of nanoelectronics and due 
to their special electronic and physical properties 
they exhibit several novel properties like high car- 
rier mobility, 3 unconventional quantum Hall effect 5 
and many others. The high carrier mobility in 
graphene demonstrates the idea for fabrication of 
high speed switching devices those have widespread 
applications in different fields. Recently GNRs are 
also used extensively in designing field-effect tran- 
sistors and this idea has been predicted in some re- 
cent nice papers. 6-8 It has many potential applica- 
tions and provides a huge interest in the community 
of nanoelectronics device research. Not only that, 
GNRs can be used to construct MOSFETs which 
perform much better than conventional Si MOS- 
FETs. In other experiment 9 it has been proposed 
that a narrow strip of graphene, the so-called a 
graphene nanoribbon, exhibits semiconducting be- 
havior due to its edge effects, unlike carbon nan- 
otubes of larger sizes which are mixtures of both 
metallic and semiconducting materials. The rea- 
son is that, in a narrow graphene sheet a band gap 
appears across the energy E = 0, while the gap 
gradually disappears with the increase of the size 
of the ribbon. It reveals a transformation from the 
semiconducting to the metallic material. This phe- 
nomenon has been studied in detail in a very recent 
theoretical work by the same author of this paper. 10 
The situation becomes quite different for the case of 
honeycomb lattice ribbons with zigzag edges. Our 
study predicts that for a ribbon with zigzag edges 
always there exists an energy gap in the conduc- 
tance spectrum across the energy E — 0. The gap 
gradually decreases with the increase of the size of 
the ribbon (as expected), but it never vanishes even 
for much larger ribbons. This clearly reveals that 
a honeycomb lattice ribbon with zigzag edges al- 
ways exhibits the semiconducting behavior, unlike 
the lattice ribbons with armchair edges where both 
the semiconducting and the metallic phases are ob- 
served by tuning the size of a ribbon. All these 
edge effects in graphene ribbons provide many key 
informations in designing nanoclcctronic devices. 

The aim of the present paper is to provide a qual- 



itative study of electron transport in honeycomb 
lattice ribbons with zigzag edges attached to two 
semi-infinite one-dimensional metallic electrodes 
(see Fig. [1]). The theoretical description of elec- 
tron transport in a bridge system has been devel- 
oped based on the pioneering work of Aviram and 
Ratncr. 11 Later, many significant experiments 12-13 
have been done in several bridge systems to un- 
derstand the basic mechanisms underlying the 
electron transport. Though in literature many 
theoretical 14-25 as well as experimental papers 12-13 
on electron transport are available, yet lot of con- 
troversies are still present between the theory and 
experiment even today. Several controlling factors 
are there which can regulate the electron transport 
in a conducting bridge significantly, and all these 
factors have to be taken into account properly to 
understand the transport mechanism. For our il- 
lustrative purposes, here we mention some of these 
issues. 

(i) The geometry of the conducting material be- 
tween the two electrodes itself is an important is- 
sue to control the electron transmission which has 
been described quite elaborately by Ernzerhof et 
al. 26 through some model calculations. 

(ii) The coupling of the bridging material to the 
electrodes significantly controls the current ampli- 
tude across any bridge system. 27 

(iii) The quantum interference effect 27-31 of elec- 
tron waves passing through different arms of the 
bridging material becomes the most significant is- 
sue. 

(iv) The dynamical fluctuation in the small-scale 
devices is another important factor which plays an 
active role and can be manifested through the mea- 
surement of shot noise, 32-33 a direct consequence 
of the quantization of charge. 

In addition to these, several other factors of the 
Hamiltonian that describe a system also provide im- 
portant effects in the determination of the current 
across a bridge system. 

Here we adopt a simple tight-binding model to 
describe the system and all the calculations are per- 
formed numerically. We narrate the conductance- 
energy and current-voltage characteristics as func- 
tions of the lengths and widths of ribbons. Our re- 
sults clearly predicts that a honeycomb lattice rib- 
bon with zigzag edges always shows the semicon- 
ducting nature irrespective of its length and width. 

The paper is arranged in this way. Following the 
introduction (Section 1), in Section 2, we present 
the model and the theoretical formulations for our 
calculations. Section 3 discusses the significant re- 
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suits, and the summary of this work is available in 
Section 4. 



2 Model and synopsis of the 
theoretical background 

Let us start with Fig. [TJ where a honeycomb lat- 
tice ribbon with zigzag edges is attached to two 
semi-infinite one-dimensional metallic electrodes, 
viz, source and drain. As the electron transport 
properties are significantly influenced by the quan- 
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Figure f : (Color online). Schematic view of a hon- 
eycomb lattice ribbon with zigzag edges attached 
to two semi-infinite onc-dimcnsional metallic elec- 
trodes, viz, source and drain. Filled circles corre- 
spond to the position of the atomic sites. 

turn interference effects, throughout the study, we 
contact the electrodes at the two extreme ends of 
nanoribbons (see Fig. [IJ to make these interference 
effects uniform. 

Calculation of the conductance: 

We use the Landauer conductance formula 34-35 
to calculate the conductance g of the ribbon. At 
very low temperature and bias voltage it (g) can be 
presented in the form, 



2e 2 



-T 



(1) 



where T gives the transmission probability of an 
electron in the ribbon which can be expressed in 
terms of the Green's function of the ribbon and its 
coupling to the two electrodes by the relation, 34-35 



T - Tr \T s G r rib T D Gr lb ] 



(2) 



where G r rih and G® ib are respectively the retarded 
and advanced Green's functions of the ribbon in- 
cluding the effects of the electrodes. The parame- 
ters Ts and T d describe the coupling of the ribbon 



to the source and drain respectively, and they can 
be defined in terms of their self-energies. For the 
complete system i.e., the ribbon, source and drain, 
the Green's function is defined as, 



G= (£-HY 



(3) 



where £ = E + iS. E is the injecting energy of the 
source electron and 5 gives an infinitesimal imagi- 
nary part to £. To Evaluate this Green's function, 
the inversion of an infinite matrix is needed since 
the complete system consists of the finite ribbon 
and the two semi-infinite electrodes. However, the 
entire system can be partitioned into sub-matrices 
corresponding to the individual sub-systems and 
the Green's function for the ribbon can be effec- 
tively written as, 

{£ - H rib - E 5 - S^r 1 (4) 



where H rib is the Hamiltonian of the ribbon which 
can be written in the tight-binding model within 
the non-interacting picture like, 



Hr 



E 

<ij> 



(4< 



(5) 



In the above Hamiltonian (H r i b ), e,'s are the site 
energies, cj (q) is the creation (annihilation) oper- 
ator of an electron at the site i and t is the nearest- 
neighbor hopping integral. A similar kind of tight- 
binding Hamiltonian is also used to describe the 
two semi-infinite one-dimensional perfect electrodes 
where the Hamiltonian is parametrized by constant 
on-site potential eg and nearest-neighbor hopping 
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and 



hD-rib9D^ D _ rib arc the self-energy opera- 



integral t - In Eq. @J Ss 

Sd 

tors due to the two electrodes, where gs and go 
correspond to the Green's functions of the source 
and drain respectively. hs- r ib and hp-rib are the 
coupling matrices and they will be non-zero only 
for the adjacent points of the ribbon, and the elec- 
trodes respectively. The matrices T$ and can 
be calculated through the expression, 
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where Y, 7 s ^ D -j and are the retarded and ad- 

vanced self-energies respectively, and they are con- 
jugate with each other. These self-energies can be 
written as 36 



S(D) 
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(7) 



where As(d) are the real parts of the self-energies 
which correspond to the shift of the energy eigenval- 
ues of the ribbon and the imaginary parts Ag^j of 
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the self-cncrgics represent the broadening of these 
energy levels. Since this broadening is much larger 
than the thermal broadening, we restrict our all cal- 
culations only at absolute zero temperature. All the 
information about the ribbon-to-electrode coupling 
are included into these two self-energies. 

Calculation of the current: 

The current passing across the ribbon can be de- 
picted as a single-electron scattering process be- 
tween the two reservoirs of charge carriers. The 
current / can be computed as a function of the ap- 
plied bias voltage V through the relation, 34 

E F +eV/2 

I(V) = j J T(E)dE (8) 

E F -eV/2 

where Ep is the equilibrium Fermi energy. Here we 
make a realistic assumption that the entire voltage 
is dropped across the ribbon-electrode interfaces, 
and it is examined that under such an assumption 
the I-V characteristics do not change their qualita- 
tive features. This assumption is based on the fact 
that, the electric field inside the ribbon especially 
for narrow ribbons seems to have a minimal effect 
on the conductance-voltage characteristics. On the 
other hand, for quite larger ribbons and high bias 
voltages the electric field inside the ribbon may play 
a more significant role depending on the internal 
structure and size of the ribbon, 36 but the effect 
becomes too small. 

3 Numerical results and dis- 
cussion 

To have a better insight in the present problem 
i.e., the dependence of the electron transport on 
the lengths and widths of ribbons, we focus our at- 
tention only on the perfect systems rather than any 
disordered one. To get these perfect systems we 
fix the on-site energies e.j = for all the sites i of 
the honeycomb lattice ribbons. The values of the 
other parameters are taken as follow. The nearest- 
neighbor hopping integral t in the ribbon is set to 
1, the on-site energy e and the hopping integral to 
for the two electrodes are fixed to and 1 respec- 
tively. The parameters t$ and Tp> are set as 0.75, 
where they correspond to the hopping strengths of 
the ribbon to the source and drain respectively. In 
addition to these, to describe the size of a ribbon 
we introduce two other parameters N and M where 



they correspond to the width and length of the rib- 
bon respectively. Thus, for example, a nanorib- 
bon with N = 1 and M = 3 represents a linear 
chain of three hexagons. Hence the parameter M 
determines the total number of hexagons in a single 
chain. Following this rule, a nanoribbon with N = 3 
and M = 4 corresponds to three linear chains at- 




Figure 2: (Color online). Conductance g as a func- 
tion of the energy E for some lattice ribbons with 
fixed width N = 3 and varying lengths where (a) 
M = 6 and (b) M = 8. 

tached side by side (see Fig. [I) where each chain 
contains four hexagons. For simplicity, throughout 
our study we set the Fermi energy Ep = and 
choose the units where c = e = h = 1. 

Let us begin our discussion with the variation of 
the conductance g as a function of the energy E. 
In Fig. [2] we plot the conductance-energy charac- 
teristics for some typical lattice ribbons with fixed 
width N = 3 and varying lengths where (a) and (b) 
correspond to the lengths M = 6 and 8 respectively. 
The sharp resonant peaks in the conductance spec- 
tra are observed almost for all energies, while for 
some other energies either the conductance g gets 
much small value or drops to zero. At the reso- 
nances the conductance gets the value 2, and there- 
fore, the transmission probability T becomes unity 
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since the relation g = 2T follows from the Lan- 
dauer conductance formula, Eq. [1] with e = h = 1. 
Now the reduction of the transmission probability 
(T < 1) for some particular energies can be ex- 
plained by considering the quantum interference ef- 
fects of the electronic waves passing through the 
different arms of the ribbon. The interpretation is 
as follow. During the motion of the electrons from 




Figure 3: (Color online). Conductance g as a func- 
tion of the energy E for some lattice ribbons with 
fixed length M = 4 and varying widths where (a) 
N = 5 and (b) N = 10. 

the source to drain through the lattice ribbon, the 
electron waves propagating along the different pos- 
sible pathways can get a phase shift among them- 
selves, according to the result of quantum interfer- 
ence. Therefore, the probability amplitude of get- 
ting an electron across the ribbon either becomes 
strengthened or weakened. This causes the trans- 
mittance cancellations and provides anti-resonances 
in the conductance spectrum. Thus it can be em- 
phasized that the electron transmission is strongly 
affected by the quantum interference effects, and 
hence the ribbon to electrode interface structure. 
Now all these resonant peaks are associated with 
the energy eigenvalues of the ribbon, and accord- 
ingly, we can say that the conductance spectrum 



manifests itself the electronic structure of the rib- 
bon. Due to the large number of energy levels, 
associated with the size of the ribbons, the reso- 
nant peaks almost overlap with each other and form 
quasi-continuous spectra in the g-E characteristics. 
The most significant feature observed from the spec- 




Figure 4: (Color online). Variation of the central 
energy gap SE as a function of the width N for some 
lattice ribbons with fixed lengths M. The red and 
blue curves correspond to M = 2 and 4 respectively. 

tra, given in Fig. [2J is that a central energy gap 
(SE) appears across the energy E = 0. With the 
increase of the length of the ribbon some more reso- 
nant peaks appear around E = 0, and accordingly, 
the gap decreases which is clearly observed from 
Figs. [21(a) and (b). Thus for a fixed width, the cen- 
tral gap can be controlled by tuning the length of 
the ribbon. 

In the same fashion, to characterize the depen- 
dence of the electron transport on the widths of the 
ribbons for a fixed length, in Fig. [3] we plot the re- 
sults for some typical lattice ribbons considering the 
length M = 4. The spectra shown in Figs. [3ja) and 
(b) correspond to the ribbons with widths N = 5 
and 10 respectively. Quite similar to the above case, 
here also a gap appears across E = and it de- 
creases with the increase of the width of the ribbon 
though the reduction is much small. Thus from the 
results described in Figs. [21 and [3] wc can empha- 
size that the width of the central energy gap always 
decreases with the increase of the size (length and 
width) of the ribbon. 

To illustrate the dependence of the gap SE with 
other system sizes, in Fig. 01 we display the varia- 
tion of SE as a function of the width N for some 
typical lattice ribbons with fixed lengths. The red 
and blue curves correspond to the lengths M = 2 
and 4 respectively. Quite interestingly we see that, 
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the gap SE gradually decreases with the increase 
of the width N, and beyond a certain value of N, 
the rate of decrease of this gap becomes much small 
and eventually it (SE) becomes almost a constant. 
Quite similar feature is also observed if we plot the 
variation of the energy gap as a function of the 
length M keeping the width N as a constant, and 
due to the obvious reason we do not plot these re- 
sults further in the present description. These re- 
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Figure 5: (Color online). Current I as a function 
of the bias voltage V for some lattice ribbons with 
fixed width N = 3 and varying lengths where (a) 
M = 3 and (b) M = 6. 

suits provide us an important signature which con- 
cern with the variation of the energy gap by tuning 
the size of the ribbon, and we can emphasize that a 
honeycomb lattice ribbon with zigzag edges always 
exhibits the semiconducting (finite energy gap) be- 
havior. 

All these basic features of electron transfer can 
be much more clearly explained from our investi- 
gation of the current-voltage (I-V) characteristics 
rather than the conductance-energy spectra. The 
current / is determined from the integration proce- 
dure of the transmission function (T) (see Eq. , 
where the function T varies exactly similar to the 
conductance spectra, differ only in magnitude by a 



factor 2, since the relation g = 2T holds from the 
Landauer conductance formula (Eq. [1}. As an il- 
lustration, in Fig. [5j we present the current-voltage 
{I-V) characteristics for some lattice ribbons with 
fixed width N = 3 and varying lengths where (a) 
and (b) correspond to the lengths M = 3 and 6 re- 
spectively. In the same footing, in Fig. [SI we plot 
the variation of the current J as a function of the 
bias voltage V for some typical lattice ribbons keep- 
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Figure 6: (Color online). Current I as a function 
of the bias voltage V for some lattice ribbons with 
fixed length M = 4 and varying widths where (a) 
N = 2 and (b) N = 3. 

ing the length as fixed {M — 4) and vary the widths, 
where (a) and (b) represent the ribbons with widths 
N = 2 and 3 respectively. The sharpness in the I- 
V characteristics and the current amplitude solely 
depend on the coupling strengths of the ribbon to 
the side attached electrodes, viz, source and drain. 
It is observed that, in the limit of weak coupling, 
defined by the condition t$(d) «t, current shows 
staircase like structure with sharp steps. While, in 
the strong coupling limit, described by the condi- 
tion Tg/m ~ t, current varies quite continuously 
with the bias voltage V and achieves large current 
amplitude compared to the weak-coupling limit. All 
these coupling effects have already been explained 
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in many theoretical as well as experimental papers 
in the literature. The key feature observed from 
these I-V characteristics is that for all such ribbons 
the electron starts to conduct beyond some finite 
bias voltage, the so-called threshold bias voltage 
Vth , and it decreases very slowly with the change of 
the size of the ribbon. Our study reveals that the 
threshold bias voltage never drops to zero, even for 
much larger systems, and accordingly, we can pre- 
dict that a honeycomb lattice ribbon with zigzag 
edges exhibits only the semiconducting behavior. 

4 Concluding remarks 

To summarize, we have addressed electron trans- 
port properties in honeycomb lattice ribbons with 
zigzag edges attached to two semi-infinite one- 
dimensional metallic electrodes within the tight- 
binding framework. We have numerically computed 
the conductance-energy and current- voltage charac- 
teristics concerning the dependence on the lengths 
and widths of the ribbons. The results have pre- 
dicted that for such ribbons i.e., ribbons with zigzag 
edges a central energy gap always exists across the 
energy E = in the conductance spectrum. The 
gap decreases gradually with the increase of the 
size of the ribbon but it never vanishes. This phe- 
nomenon clearly manifests that a honeycomb lat- 
tice ribbon with zigzag edges always shows semi- 
conducting nature, unlike the lattice ribbons with 
armchair edges where both the semiconducting and 
the metallic phases are observed by controlling the 
size of the ribbon. This semiconducting behavior of 
the lattice ribbons with zigzag edges has been much 
more clearly addressed from our presented current- 
voltage characteristics. It has been observed that 
the current starts to appear beyond some finite bias 
voltage i.e., the threshold bias voltage V t h has a 
non-zero value. This Vth doesn't change apprecia- 
bly with the increase of the size of the ribbon and we 
have also examined that the threshold bias voltage 
never reduces to zero even for much larger systems, 
which predicts the semiconducting nature only. 

This is our first step to describe how the elec- 
tron transport properties in honeycomb lattice rib- 
bons with zigzag edges depends on the size of the 
ribbons. We have made several realistic assump- 
tions by ignoring the effects of the electron-electron 
correlation, disorder, interaction with a substrate, 
temperature, finite width of the electrodes, bound- 
ary of the ribbons, etc. Here we discuss very briefly 
about these approximations. The inclusion of the 
electron-electron correlation in the present model 



is a major challenge to us, since over the last few 
years people have studied a lot to incorporate this 
effect, but no such proper theory has yet been devel- 
oped. In this work, we have presented all the results 
only for the ordered systems. But in real samples, 
the presence of impurities will affect the electronic 
structure and hence the transport properties. Be- 
side these, in experiments, the graphene nanoribbon 
is deposited on an insulating substrate which has 
also not been included in our present study, and, 
it has been observed from first-principles calcula- 
tions that the effect of the substrate is too week. 37 
The effect of the temperature has already been 
pointed out earlier, and, it has been examined that 
the presented results will not change significantly 
even at finite temperature, since the broadening of 
the energy levels of the ribbon due to its coupling 
with the electrodes will be much larger than that 
of the thermal broadening. 34 The other important 
assumption is that here we have chosen the linear 
chains instead of wider leads, since we are mainly 
interested about the basic physics of the ribbon. 
Though the results presented here change with the 
increase of the thickness of the leads, but all the ba- 
sic features remain quite invariant. The effect of the 
boundary is also an important issue in this context. 
Here we have considered only the perfect geometry 
of the nannoribbons. Several interesting features 
will be observed for the nanoconstrictions with dif- 
ferent shapes 38 like, square-shaped, wedge-shaped 
nanoconstrictions, etc. Finally, we would like to 
say that we need further study in such systems by 
incorporating all these effects. 
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